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We measure the top quark mass (mt) in pp collisions at a center of mass energy y/s = 1.96 TeV 
using dilepton — >■ W^hW~h — >■ £'^V(b£~ i>eb events, where £ denotes an electron, a muon, or a tan 
that decays leptonically. The data correspond to an integrated luminosity of 5.4 fb~^ collected with 
the DO detector at the Fermilab Tevatron Collider. We obtain rut = 174.0±1.8(stat)±2.4(syst) GeV, 
which is in agreement with the current world average nit = 173.3 ± 1.1 GeV. This is currently the 
most precise measurement of mt in the dilepton channel. 

PACS numbers: 14.65. Ha 



The measurement of the properties of the top quark 
has been a major goal of the Fermilab Tevatron Col- 
hder experiments since its discovery in 1995 [H, Q- As 
the heaviest known elementary particle, the top quark 
may play a special role in the mechanism of electroweak 
symmetry breaking. A precise measurement of its mass 
(mt) is of particular importance, since, combined with 
the measurement of the W boson mass, it provides an 
indirect constraint on the mass of the Higgs boson in the 
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standard model (SM), and can also constrain possible 
extensions of the SM. 

We present a new measurement of the top quark mass 
in the dilepton channel (ee, e/z, /z/i) in tt W^hW^b 
£'^i>ib£~Dib events, where £ denotes an electron, a muon 
or a tau decaying leptonically, using the matrix element 
method. The first measurement of mt based on this 
method was performed in the lepton+jets channel by 
the DO experiment The CDF Collaboration has ap- 
plied the matrix element approach to determine mt in 
the dilepton and all-hadronic final states 0, , obtaining 
a mass precision of 4.0 GeV for dilepton events The 
measurement of mt in the dilepton channel has also been 
carried out by using other techniques reaching a 

precision of 3.7 GeV. We report a measurement based 
on data collected by the DO detector, corresponding to 
5.4 fb~^ of integrated luminosity from pp collisions at 
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= 1.96 TeV. 

The DO detector has a central tracking system, con- 
sisting of a sihcon microstrip tracker and a central fiber 
tracker, both located within a 1.9 T superconducting 
solenoidal magnet [l^ , with the design providing track- 



ing and vertexing at pseudorapidities jryj < 3 [13|. The 
liquid-argon and uranium calorimeter has a central sec- 
tion covering pseudorapidities \r]\ up to « 1.1 and two 
end calorimeters that extend coverage to \r]\ w 4.2, with 
all three housed in separate cryostats 1J|. A muon sys- 
tem outside the calorimeters covers I??] < 2 and consists 
of a layer of tracking detectors and scintillation trigger 
counters in front of 1.8 T toroids, followed by two similar 
layers after the toroids 15 1. 

Despite the small branching fraction of this final state 
and the presence of two neutrinos in each event, the mea- 
surement of TOt in the dilepton channel is interesting be- 
cause the lower background and the smaller jet multiplic- 
ity relative to the lepton-|-jets channel result in a reduced 
sensitivity to the ambiguity from combining jets in the re- 
construction of rrit ■ The dilepton measurement therefore 
complements the results from other final states. More- 
over, significant differences in measured values of rrit in 
different tt decay channels can be indicative of the pres- 
ence of physics beyond the SM [l^ • 

As the SM predicts top quarks to decay almost 100% 
of the time into a W boson and a b quark, tt events 
are classified according to the decays of the W boson. 
In the dilepton channel, both W bosons decay lepton- 
ically, — > t^fg [17[ with £ — e, fj, or t. We ana- 
lyze the events characterized by two leptons ee, e/x, or 
fifi, with a large transverse momenta (pt), large im- 
balance in transverse momentum from the undetected 
neutrinos (^t), and two high-py jets from the b quarks. 
The W'^ t^Vt- decays contribute through secondary 
— > l^vgD-r transitions. For the ee and ^/z analysis, 
we consider events selected by a set of single-lepton trig- 
gers. For the e/i channel, we use a mixture of single 
and multilepton triggers and lepton-|-jet triggers. Dilep- 
ton tt events are required to have at least two oppositely 
charged, isolated leptons with pT > 15 GeV, and either 
|ry| < 1.1 or 1.5 < \ri\ < 2.5 for electrons and \r]\ < 2 
for muons. If more than one lepton-pair combination is 
found in an event, only the pair with the largest sum 
in scalar pt is used. Events must have at least two 
jets with Pt > 20 GeV and \r]\ < 2.5, well separated 
from the selected electrons. No explicit 6-jet identifica- 
tion is required in this analysis. The main sources of 
background in the dilepton channel are Drell-Yan and 
Z boson production {Zj^* — > I'^l^), diboson produc- 
tion {WW, WZ, ZZ), and instrumental background that 
originates from limited detector resolution and lepton 
misidentification. In the ee channel, the discrimination 
between the ti signal and background improves by re- 
quiring a large significance of the measured jix, which 
is defined through a likelihood discriminant constructed 



from the ratio of to its uncertainty 18|. In the niJ, 
channel, we require, in addition, > 40 GeV. In the 
efi channel, the requirement Ht > 115 GeV, where Ht 
is the scalar sum of the transverse momenta of the lead- 
ing lepton and the two leading jets, rejects most of the 
contribution from t+ — > t'^viVr- The above selections 
minimize the expected statistical uncertainty on toj. In 
total, we select 479 candidate events with 73, 266, and 
140 events, respectively, in the ee, e/i, and ///i channels, 
of which about 13 ± 5, 48 ± 15, and 56 ± 15 events, re- 
spectively, are expected to arise from the background. 

The matrix element method is based on the probabil- 
ity for a given event to resemble a signal, which depends 
on the value of toj, or a background, which is usually 
independent of nit. Assuming that the different physics 
processes leading to the same final state do not inter- 
fere, the event probability can be written as the sum of 
probabilities from all possible contributions. In practice, 
because the matrix element method requires significant 
computing time, only the dominant background is taken 
into account, and the total event probability is given by 



htPttix; mt) + (1 - fu)Pz+2 jcts(a;), (1) 



where fti is the fraction of tt events, Pu and Pz+2 jets 
are the signal and background probability densities, re- 
spectively, mt is the assumed top quark mass, and x 
reflects the observed kinematic variables, i.e., the four- 
momenta of the measured jets and leptons. In the ee, 
/i//, and e^ channels, Z + 2 jets events with Z — >■ e"'"e~, 
Z — > ^jL'^fi~ and Z — > r+r^ — >■ e'^Ue^jL~iy^ are the dom- 
inant source of background. The second leading back- 
ground, from misidentified leptons, is approximately a 
factor of 3 smaller. While neglecting the other back- 
ground probabilities leads to some bias, the calibration 
procedure described below allows us to correct for these 
and other limitations of the model. 

The leading-order (LO) matrix element for qq ^ tt ^ 
W~^bW~b — >■ i'^vib£~V(b is used to compute the tt prob- 
ability density. For each final state y of the six produced 
partons, the signal probability is given by 



Ptt(x;mt) 



1 



Cobs ("It) 



f 

^ / dqi dq2 /pDF(gi) /pDFfe)- 
■i—i 



(2) 



(2^)4|M(y;™0l' 



d$6 W{x,y) W{p*S) , 



where qi and q2 denote the momentum fractions of 
the incident quarks in the proton and antiproton, respec- 
tively, /pDF are the parton distribution functions (PDF) 
for finding a parton of a given flavor and longitudinal mo- 
mentum fraction in the proton or antiproton (in this anal- 
ysis we use the CTEQ6L1 PDF 0), s is the square of 
the energy in the qq rest frame, M{y) is the leading-order 
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matrix element [2^ and d$6 is an element of the 6-body 
phase space. Detector resolution is taken into account 
through a transfer function W{x,y) that describes the 
probability of the partonic final state y to be measured 
as X. The finite transverse momentum of the tt system 
is accounted for through an integration over its probabil- 
ity distribution, which is derived from parton-level sim- 
ulated events using ALPGEN [21|, employing pythia [22 | 
for parton showers and hadronization. As the angular 
resolution of the jets and leptons, as well as the elec- 
tron energy resolution, are sufiiciently well determined, 
there is no need to introduce resolution functions. By 
taking into account energy and momentum conservation, 
Eq. ([2]) can be reduced to an integration over the en- 
ergies associated with the b quarks, the lepton-neutrino 
invariant masses squared, the differences between neu- 
trino transverse momenta, the transverse momentum of 
the tt system, and the radii of curvature (p^^) of muons. 
The sum runs over both possible jet-parton assignments 
and over up to two real solutions for each neutrino en- 
The normalization factor Cobs is the product 
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ergy 

of the LO cross section and the mean efficiency of the 
final selections. A transfer function W{x, y) is used for 
each jet and each muon in the final state. The jet energy 
resolution is parametrized as the sum of two Gaussian 
functions, with parameters depending linearly on parton 
energies, while the resolution in muon is described 
by a single Gaussian function. All parameters in W{x, y) 
are determined from Monte Carlo (MC) tt events, tuned 
to match the resolutions observed in the data. 

To take account of all background processes and to 
provide a correct statistical sampling of possible spin, fla- 
vor, and color configurations, the background probability 
Pz+2 jets is calculated by using VECBOS ^ . Since Z 
T^T~ decay is not modeled in vecbos, an additional 
transfer function in the efi channel is used to describe 
the energy of the final state lepton relative to the initial 
T lepton, derived from parton-level information (23| . The 
direction of the final state lepton is assumed to be close 
to that of the r lepton, since only in such cases is the 
lepton from the r decay sufiiciently energetic to pass the 
Pt selection. For the {Z — > t+t^ — > e+i^e/i^i/p) -I- 2 jets 
probability, the energy fractions for final state leptons 
are sampled according to this r transfer function. The 
jet and charged-lepton directions are assumed to be well- 
measured, and each kinematic solution is weighted ac- 
cording to the Pt of the Z + 2 jets system. The integra- 
tion of the probability for Z + 2 jets is performed over the 
energies of the two partons that lead to the jets. Both 
possible assignments of jets to quarks are considered. 

To calculate the signal and background probability 
densities, a MC-based integration of Eq. ([2]) is performed 
and rrit is changed in steps of 2.5 GeV over a range of 
30 GeV. For each mass hypothesis, a likelihood func- 
tion Ltot {fTT-t , ftt) is defined by the product of individual 
event probabilities Pcvt , and the signal fraction f is de- 



termined by minimizing — In Ltot • Finally, the most likely 
value of rrit and its uncertainty are extracted from a fit of 
Ltot{mt) to a Gaussian form near its maximum by using 
the value of fu found in the previous step. 

To check for any bias caused by approximations of the 
method, such as the use of the LO matrix element for 
Ptt or from neglecting backgrounds other than Z -I- 2 jets, 
the measurement is calibrated by using MC events gen- 
erated with ALPGEN -I- PYTHIA. All evcnts are processed 
through a full GEANtS [25| detector simulation, followed 
by the same reconstruction and analysis chain as used for 
the data. Effects from additional pp interactions are sim- 
ulated by overlaying the data from random pp crossings 
over the MC events. Five ti MC samples are generated 
with input top quark masses of mt = 165, 170, 172.5, 
175, and 180 GeV. Probabilities for the ti signal and for 
Z/7* — 7> diboson and instrumental backgrounds, 

are used to form randomly drawn pseudoexperiments. 
The total number of events in each pseudoexperiment is 
fixed to the number of events in the data for the combined 
dilepton channels. The signal and background fractions 
are fiuctuated according to multinomial statistics around 
the fractions determined from the measured ti cross sec- 
tion in the separate channels [2^. The mean values of 
mt measured in 1000 pseudoexperiments as a function of 
the input mt are shown in Fig. [Ija). The deviation from 
the ideal response, where the extracted mass is equal to 
the input mass, is caused both by the presence of back- 
grounds without a corresponding matrix element in the 
event probability and by approximations in the calcu- 
lation of the Z + 2 jets probabilities. For the case of 
background-free pseudoexperiments, no difference is ob- 
served. The width of the distribution of the pulls ("pull 
width"), defined as the mean deviation of mt in single 
pseudoexperiments from the mean for all 1000 values at 
a given input mt, in units of the measured uncertainty 
per pseudoexperiment, is shown in Fig. [Ijb). The sta- 
tistical uncertainty measured in the data is corrected for 
the deviation of the pull width from unity. The calibrated 
value of mt from the fit to the data is shown in Fig.[2][a). 
Figure EJ^b) compares the measured uncertainty for mt 
with the distribution of expected uncertainties in pseu- 
doexperiments at mt = 175 GeV. The difference between 
the observed and median expected uncertainty is not sta- 
tistically significant. We also note that, when we change 
the signal to background ratio within uncertainties, the 
expected uncertainty generally increases and agrees well 
with the observation. 

Systematic uncertainties on the measurement of mt 
can be divided into three categories. The first involves 
uncertainties from modeling of the detector, such as the 
uncertainty on the energy scale of light-quark jets and 
the uncertainty in the relative calorimeter response to h 
and light-quark jets, as well as in the energy resolution 
for jets, muons, and electrons. The second category is re- 
lated to the modeling of production. This includes pos- 
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sible differences in the amount of initial and final state ra- 
diation, effects from next-to-leading-order contributions 
and different hadronization models, color reconnection, 
and modeling of 6-quark fragmentation as well as uncer- 
tainties from the choice of PDF. The third category com- 
prises effects from calibration, such as the uncertainties 
in the calibration function shown in Fig. [Ija), and from 
variations in signal and background contributions in the 
pseudoexperiments. Contributions to the total system- 
atic uncertainty in the measurement of nit are summa- 
rized in Table HI 

The dominant systematic uncertainty arises from the 
different detector response of light and 6-quark jets. It 
accounts for the different calorimeter response of single 
pious in the data and MC simulation and the different 
fractions of single pions in light and &-quark jets. The 
relative uncertainty of the response has been evaluated 
to be 1.8% leading to a shift of 1.6 GeV in m^. The next 
important uncertainty comes from uncertainties in the jet 
energy scale (JES) of light quarks. This JES is calibrated 
by using 7-l-jets events More than 80% of the JES 
uncertainty is due to the understanding of the detector 
response and the showering of jets. The total uncertainty 
typically adds up to about 1.5% per jet, which trans- 
lates into an uncertainty on rrit of 1.5 GeV. The main 
uncertainty from modeling tt production is from higher- 
order effects and hadronization. It is evaluated by using 
tt events generated with MC@NLO [l^l and evolved in 
HERWIG [29|. The next leading uncertainty on modeling 
tt arises from the description of 6-quark fragmentation. 
It is derived by comparing the extracted mt for the de- 
fault measurement with the result using a reweighting 
of the default MC samples to a Bowler scheme tuned to 
LEP or SLD data fsoj. The largest difference is quoted 
as the uncertainty. 

In summary, we have presented a measurement of the 
top quark mass in the tt — >■ W~^bW~b (,^viM^V(h 
channel using the matrix element method. Based on an 
integrated luminosity of 5.4 fb^^ collected by the DO Col- 
laboration, the top quark mass is found to be 



mt = 174.0 ± 1.8(stat) ± 2.4(syst) GeV. 



(3) 



This measurement is in good agreement with the cur- 
rent world average mt — 173.3 ± 1.1 GeV (sij. Its total 
uncertainty of 3.1 GeV corresponds to a 1.8 % accuracy 
and represents the most precise measurement of mt from 
dilepton tt final states. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST (In- 
dia); Colciencias (Colombia); CONACyT (Mexico); KRF 
and KOSEF (Korea); CONICET and UBACyT (Ar- 
gentina); FOM (The Netherlands); STFC and the Royal 
Society (United Kingdom); MSMT and GACR (Czech 



> 

O 



3 -5 

CO 
0) 







/ 


uu 




/ 

/ 


(a) 




/ . 

/ / 




/ r 




A 

//^ 
/y 




/ 
/J 
// 
// 
/ K 
/ , 


/ X 




' offset 


: -1.36 ±0.10 


slope 


+0.97 ± 0.02 








-5 5 

Input m^-172.5 (GeV) 



-5 5 

Input m -172.5 (GeV) 



FIG. 1: (a) Mean values of mt and (b) pull width from sets 
of 1000 pseudoexperiments as a function of input mt for the 
combined dilepton channels. The dashed lines represent the 
ideal response in (a), where the extracted mass is identical to 
the input mass, and in (b), where the statistical uncertainty 
requires no correction. 
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FIG. 2: Combined for all channels: (a) Calibrated and nor- 
malized likelihood for the data as a function of mt with 
the best estimate as well as 68% confidence level region 
marked by the shaded area and in (b) the expected distri- 
bution of uncertainties with the measured uncertainty indi- 
cated by the arrow. As the top quark mass is measured to 
be mt = 174.0 GeV, the expected distribution in (b) is shown 
for the closest input mass mt = 175 GeV used in the pseudo- 
experiment. 
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